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ON THE PARTICULAR INTEGRALS OF THE PRANDTL-BUSEMANN ITERATION EQUATIONS 

FOR THE FLOW OF A COMPRESSIBLE FLUID ‘ 

By Carl Kaplan 


SUMMARY 

The ‘particular integrals of the second-order and third-order 
Prandtl-Busemann iteration equations for the flow of a com- 
pressible fluid are obtained by means of the method in which 
the complex conjugate variables z and z are 'utilized as the 
independent variables of the analysis. The assumption is made 
that the Prandtl-Glauert solution of the linearized or first-order 
iteration equation for the two-dimensional flow of a compressible 
fluid is knoum. The forms of the particular integrals, derived 
for subsonic flow, are readily adapted to supersonic flows with 
only a change in sign of one of the parameters of the problem. 


systematic manner. The eipKcit expressions obtained for 
the particular integrals of the second- and third-order itera- 
tion equations are believed to yield essentially the solution 
of the problem of high subsonic flow past an arbitrary two- 
dimensional profile, since it is never a difficult problem to 
supply the solutions of the homogeneous equation necessary 
for the fuiflUment of the boundary conditions. It is notet- 
worthy that the particular int^rals, derived for subsonic 
flow, can be adapted to supersonic flow with simply a change 
in sign of one of the parameters. 

FUNDAMENTAL EQUATIONS 


INTRODUCTION 

For the past several years iteration methods have been 
increasingly applied to the solution of compressible-flow 
problems. The most useful method from the point of view 
of aeronautical applications and the one discussed in this 
report is based on small perturbations with respect to the 
undisturbed flow. The Prandtl-Glauert and Ackeret solu- 
tions in two-dimensional subsonic and supersonic flow, 
respectively, obtained by means of the linearization of the 
fundamental nonlinear differential equation for cdmpressible 
flow, are presumed to be known and are taken as the initial 
steps in this iteration process. Higher-order solutions are 
then obtained by retaining appropriate powers and products 
of the perturbation quantities. This method of iteration 
has been variously labeled the Ackeret iteration process and 
the Prandtl-Busemann small perturbation method when 
limited to two-dimensional subsonic flow. The procedure 
has been extended in recent years to both two-dimensional 
and axisymmetrical supersonic-flow problems. 

In a recent publication (reference 1), Van Dyke succeeded 
in obtaining by trial the particular integral of the nonhomo- 
geneous second-order iteration equation for the velocity 
potential in supersonic flow. The general solution is then 
easily obtained by adding solutions of the homogeneous 
equation with proper r^ard to the boundary conditions at 
the surface of the solid and at infinity. 

The purpose of the present report is to show a procedure 
by means of which the partictilar integrals of the higher- 
than-first-order iteration equations can be derived in a 


The fundamental nonlinear differential equation governing 
the flow of a compressible fluid is 

(&+!?)=“ 

where 

X, T rectangular Cartesian coordinates in flow plane 

u, V fluid velocity components along X- and V-axis, respec- 
tively 

c local speed of sound 

The condition for irrotational motion is that _ 

and leads to a velocity potential # defined by 


'^~bX 



( 2 ) 


If the body is held fixed in a uniform stream of velocity U, 
the relation between the local speed of sound c and the speed 
of the fluid is given for adiabatic processes by 



Mj(^l 


u^+v^\ 
IP ) 


( 3 ) 


> St^)ers(;des NAG A T^^ 2159, “On the Partfcular Integrals of the Prandtl-Bosamann Iteration EquaBcms Cor the Flow of a Compressibla Flold" hy Cari Kaplan, 1950. 
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where 


sound speed in undisturbed fluid 

y ratio of specific heats at constant pressure and constant 
volume 


Mach number of undisturbed stream (Ujc„) 


With the introduction of a characteristic length I as unit 
of length and the undistiubed stream velocity U as unit of 
velocity, the quantities X, Y, u, v, and $ for the remainder 
of the analysis denote, respectively, the nondimensional 
quantities X/l, Yjl, ujU, v/U, and while c and retain 
their original meanings. By means of equations (2) , equa- 
tions (1) and (3) then become, respectively, 

^xz+(j-i—MJv^^^rr—2MJ^z^y^xr=0 

(4) 

and 


£i= 1 MJ[1 -(uH »*)] (5) 


where the subscripts X and Y denote partial differentiations 
with respect to the designated variables. 

In order to obtain the iteration equations based on small 
perturbations of the imdistuxbed stream, the assumption is 
made that the velocity potential $ can be expanded in the 
form 

... ( 6 ) 

For the purpose of defining and controlling the iteration 
procedure, the function and its derivatives are then 
regarded as small compared with the preceding approxima- 
tion and its derivatives. 

From equations (2) and (6), 


■W=^.r=l + '^’ix + ^2X+^&r+ • • • 
and 

When these expressions for u, v are introduced into equation 

(4) , together with the expression for cVc»” given by equation 

(5) , and the powers and products of €>» and their derivatives 
are grouped according to the assumptions of the small per- 
turbation method, the following iteration equations for the 
first three approximations $i, and result: 


(1 — M 1^2 [(7+ iK’ixxd' 

(7 — l)$iyy]-l-2 + 

[(7 + l)^sxx + (7 ~ 1 )^2y y] + 1 y^ixy + 

^jy^ixy -(-$iy$12ty)| (9) 

For slender bodies, the first few steps of this iteration process 
may be expected to yield an accurate result with the exceji- 
tion of a small region in the neighborhood of a stagnation 
point. Even at stagnation points, the iteration method has 
been shown to represent correctly the effect of compressibility 
(references) . The accuracy of the calculations obviously de- 
pends upon the number of terms determined, each additional 
term reducing the regipn of inaccuracy in the neighborhood 
of a stagnation point. 

The iteration equations (7), (8), and (9) may be put into 
more familiar forms by the introduction of a new set of 
independent variables x and y, where 

a:=A' ) 

, > flO) 

l/=FVr^) 

Thus, for M„<1, equation (7) is transformed into a Laplace 
equation; whereas equations (8) and (9) are transformed into 
Poisson equations with the right-hand sides composed of, 
respectively, double products and triple products of previ- 
ously determined perturbation quantities. It is further as- 
sumed that the solution of equation (7) is available. This 
initial step in the approximation to the e.xact nonlinear solu- 
tion is usually easily obtained, as it represents the Prandtl- 
Glauert approximation (reference 3, appendix B). The 
purpose of the present report is then to derive explicit ex- 
pressions for the particular solutions of the second- and 
third-order iteration equations (8) and (9). 

CALCULATION OF THE PARTICULAR INTEGRAL OF THE 
SECOND-ORDER ITERATION EQUATION 

By introducing the independent variables x and y defined 
by equation (10), the second-order iteration equation (8) 
becomes 


where 




7+1 MJ 
2 |8* 


( 11 ) 


(1 — M„^)$ixx+'^iyr = 0 (7) 

(1 — MJ)^iXX +^2Y Y = MJ[{y+ l)'i'ix$ixx + 

(7 — l)$ix'l’iyy + 2 #iy$ix-y] ( 8 ) 


and where use has been made of Laplaee’s equation 


'^1x1+ '^’iw — 0 
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The procedure for obtaining the particular integrals of the 
higher-order iteration equations is based on the use of the 
complex conjugate variables 2 and ~z as independent variables. 
Thus, 

z = x-\riy 
z=x—iy 

and the equivalence of operators 

A=A+A 

Si dz~dz 


A— • / ^ ^ \ 

Qy ^ \d 2 dzj 


a* a* 
ar® as* 


.2^-lA 

^ bzbz^b-p 


^ a* „ a® a® 

ay® as® as as as® 

Then equation (7) for becomes 

The most general real solution of this equation is 

[Wi(s)-i-¥i( 2 )l (13) 

or 

€>i = R.P. Wi(z) = R . P. Wi(z) 

where tOi(s) is an arbitrary analytic fimction of s, tOi(s) is its 
conjugate complex, and where the symbol R.P. stands for 
“real part of.” The imaginary part of a’i(s) is a fvmction fi, 
say, related to #1 by means of the Cauchy-Riemann equations 
and hence also satisfies Laplace’s equation. The function 
does not represent the stream function of the actual com- 
pressible fiow and does not appear in the final expressions 
of the particular integrals. The following relations will be 
found useful and are easily verified: 


$ 1=2 («?i+Wi) 


^lr=2 (Wu + Wit) 


®lir=o (Wli—Wld 


^lir=7^ — Will) 


<&Irt=o (Wi„d-Wiif) 


Then 

2$u'l'to=R- P- (Wi, + iCiy) Wiu 

—R.P. (V}u—Wu)Wuz 

■ and equation ( 11 ) for the second approximation #4 becomes 
MJR.'P. [<7WuWi„-l-(2-i-o-)iFiTWi«] 

=i MJ R. P. [1 (wu*),+(2 -H v) (wuTcii),] 

If $2 is defined to be the real part of a nonanalytic function 
Wi(z,z), then 

MJ [f (wu*),-K2 -i- v) (wi.wu),] (14) 

This equation can be integrated immediately by inspection 
and yields the general solution 

MJ [1 zwJ+(2 + a) WiWu+F(z)j (16) 

where, because only the real part $2 of its is of interest, only 
one arbitrary analytic fvmction F(s) need be included. The 
function F(z) satisfies Laplace’s equation and is so chosen 
as to satisfy the required boundary conditions at the surface 
of the body and at infinity. The part of the expression on 
the right-hand side of equation (15), excluding the arbitrary 
function F(s), is the particular integral of equation (14) and 
may be expressed in real form in the following manner: 
Suppose 

F(s)=— I sWi,*-i-(2-l-(r)wii0u-b/(z) (16) 

where /(s) is again an arbitrary analytic function of s. 

Then with the aid of the relation 

where use has been made of the Cauchy-Riemann conditions 

^ir=—4'ix 

the expression for 4>i, obtained from equation (15), becomes 

^r=\MJ [- 2 < 7 y$„,$i,-f 2(2-H(r)$i#i,4-R.P./(2)] (17) 

The expression on the right-hand side of this equation, 
excluding R.P./(s), namely, 

^,=MJ [(l -hf) -fi-i «i>ix 


(18) 



912 


REPORT 1030 — NATIONAL ADVISORY COMMITTEB FOB AERONAUTICS 


corresponds precisely to the particular integral obtained by Van Dyke (reference 1) for two-dimensional supersonic flow 
with ff replaced by — o-, where for supersonic flow the deflnition of a is 

r+1 MJ 

2 . ■ 

It is rather noteworthy that the particular integral of the second-order iteration equation (8) can be obtained for both sub- 
sonic and supersonic flows by simply interchanging the sign of the parameter a. 

CALCULATION OF THE PARTICULAR INTEGRAL OF THE THIRD-ORDER ITERATION EQUATION 

In this section, the particular integral of equation (9) involving only $i, $2, and their derivatives is derived. For this 
purpose, the variables x,y and the parameter <r are introduced. Equation (9) then takes the following form: 

Use is again made of the complex conjugate variables z and i as independent variables. Thus, 

^i=R.F.Wi(z,'F) 

<i>2*=R. P. (-wiij-f-toai) 

<i>2,=R . P. i (i02z — W25) 

^2zx~ R . P . (Wa ,, 2 Wm -|- W az) 

#2,»=R. P. i(Wu,—Wat). 

P. (Wu-l-Wij)‘Wi„ 

-R. P- J (Wii- 

1 L — 

^i«’J> 2 *+'Six<E> 2«=R. p. 2 (Wi„+Wm)(wj.+MJ2s)+R.P. jXwi>+y’ij)i'U>s»+2Wbj+Wssi) 

R.P. — WihXW 2 .— Wh) . 

—R. P. i (wu—Wa)(Wi„—W 2 n) 

— R. P. (Wi,*— Wn*)Wi„ 

Then equation (19), with#a=R-P- Wi{z,z), can be written as foUows: 

l^3«r=j(l— |9*)v(Wa*Wt„+i02*z«'i»+WjsWin+WznWij)-Hj(l— )3*)(2-h(r)(WtoWin+W2,,Wii-hWaWi»+WsnWii) + 
j^(l-/3^<7[)3*(l + 2v)-(l-2v)]wi.*tOi„+j^(l-j8*)D3^4-f5<7-h2<r*)-^a(3-f-2£r)]ii;ir*Wi„+ 
i (1 -]8“)[i3*(2 + 5 ff -h 2 (-2 + (T 2 

Introducing the expressions for vh{z,z) and its derivatives with respect to 2 and 2 from equation (15) yields for ttij,! the 
following equation: 

Wte!=^(l — lS*)*(7=^(«Ji,*Wi„)*-l-j^(l -iS7'o-(2-f-o-)(2Wi)i(wi,Wi„),-l-j^(l -i8*)*<r(wi,F0,-|-j^ (1— /3*)o-[/3®(2-l-(r)+^ 

(4-f-3v)lijis(u>i,*).-l-^ (1 — /3*)V(2-|-v)(wi,*).(2WiT)s-f ^ (1— i3*)*(2-hff)*Wi„(iiJi«)ii),-i--j^ (1 — iS*)*(2-f-ff)(wjiTF,)xi+ 

^ (1 -^*)*(2-|-<7)^(iiji Vix«+^ (I-|3*)[^*(8-hl0<r+3<r*)-l- <r(6 -(- 5 <r)]«J.„i/li.* 


( 20 ) 
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In the derivation of this egDression for Ws.-, free use was made of the fact that insofar as the real part of lo* is con- 
cerned, terms of the nature g(z)h{z) and g(z)h{z) are equivalent. It is important to note that this type of operation leaves 
unaltered the real part of Wj. Since is the quantity sought in the calculations, changes in the imaginary part ojf tOg 
are of no consequence in the final results. 

Equation (20) can be integrated immediately by inspection and yields the following result: 

j3y(72Wul^x+j^(l— j3^ff[j3*(2 + tr)-K4-|-3(r)]iFiWi.*4- 

^ (1 -^*)V(2 + <r)zw,tWu*+~ (1 -m + «r)*i0iMu+^(l-|37(2+ <r)(Ew:,-lF,iCi)-|-^(l-j8*)X2-l <y)*W,*Wr„+ _ 

^(l-/S0<ri(3/3=-t-5)2Wi,*-l-^ (l-/3*)[j8‘(8-i-10<r+3O+46 + 5<7)Iii7iiJii;u*£f2 (21) 

Equation (21) is the particular integral of equation (20). The most general solution is obtained by adding an arbitrary 
analytic function G(z), satisfying the homogeneous or Laplace equation G,t=0. An arbitrary' function of z, customarily in- 

eluded in the general solution, need not be consider^ here because only the real part of itj is of interest. In fact,^ t]^ ^ 

omitted arbitrary function is the complex conjugate G(z). 

In order to obtain the desired form of (the real part of Wt) from equation (21), F(z) is replaced by its expression ~~ 
given in equation (16), and the real part oif(z) in equation (17) is replaced by 

R. P. /(2)=-^j#aH-2 2(2 -1- 

The final form of the particular integral of the third-order iteration equation (9) then becomes 

(1 -/3*X2 + <rX#n^a+#i$ix)-l- 

|(1 -i3*)*(r*y’$,#i,$^-l-^(l -^*)^[(6 + 5^)+3(-2-f'«r),3*]j/^,,’+ 

^(1 -^*)<r[(10-<r)-(10-t-7,r)^*]y$u^i,+i(l-^^M2 + <r)|/#i($ix#i^+$„$:«)-f 

^(l-/3*)[(-16-10v-lO-Kl6-l-22<r-l-7(r^)/S^$i#i,*-^(l-|3^o-[(6-l-5<r)-l-3(2-t-<r))3^$,V- - 

|(l-)3=)*(2-f v)*$i=$ux+^(l-|3^k(6-h5v)-l-(8-l-10<r-h3<ri)|S*]<^i,J dy] (22) 

The corresponding expression for for supersonic flow is obtained by simply replacing o- by — o- with 

7+1 MJ 
2 MJ-1 

and jS® by — /3* with j3-=A/„*— 1. The physical plane variables X and Y are easily inserted into both equations (18) and 
(22) by means of the transformation equations (10), 

x=X 

y=PY 

It is pointed out that the forms of the two particular integrals, equations (18) and (22), derived in this report are identica l 
for both subsonic and supersonic flow. The apparent differences are caused by a change in sign of the parameter j8*. Thus, 
jS* and <r. are positive for both subsonic and supersonic flow. Actually, of course, the fimetions represented by $j, 

'J'j, . . . are different for the two types of flow. For subsonic flow, these functions are derived from analytic and nonanalytic 
functions of z and 2 ; whereas for supersonic flow, they involve the real “characteristics” variables z±j3y. 

Note that the last term of the e.xpression on the right-hand side of equation (22) contains the indefinite integral 

7= J [($tE=— #„®)da:+ 2 c&i^$j, dy\ 


(23) 
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It is obvious from the corresponding complex integral in 
equation (21) that the integrand of equation (23) is an exact 
differential. This fact can also be easily verified with the 
help of Laplace’s equation 


Thus, 




dy 




(24) 


Equation (24) represents the necessary and sufficient condi- 
tion that the integrand of equation (23) -be an exact differ- 
ential. Further, according to the theory of exact differential 
equations, the integral I may be expressed as follows: 


I=j^Mdz+j dy 

where 


(25) 


and where by J" M dx is meant the result of integrating 

Mdx with y considered constant. The expression within 
the parentheses, namely 


iV~ 

is a function of y only.' This statement can be verified as 


follows: Thus, 

and because y is considered constant in the process of inte- 
grating Mdx, it is clear that 


i>xj 


Mdx=M 


Hence the right-hand side of equation (26) is 


pN_ 

6x 


pM 
dy ' 


which vanishes because of the condition for the existence of 
an exact differential. 

Note that in general it is simpler to perform the complex 
integration J Wudz rather than transform to a real integral 
and then perform the integration. 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., May S9, 1960. 
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